基于混合遗传算法的多冲量最优变轨

针对航天器多冲量最优变轨问题,建立了多冲量最优变轨的数值优化模型,给出了一种遗传算法(GA)与序列二次规划算法(SQP)结合的混合优化算法.该算法不需初值猜测,全局和局部搜索能力强并且计算效率高.仿真计算了燃料最优交会和时间最短拦截问题,比较了GA,SQP以及混合遗传算法的性能.针对混合遗传算法得到的不同冲量次数变轨的优化结果,分析了冲量次数对变轨性能指标的影响.结果表明混合遗传算法综合性能最高,冲量次数对不同性能指标的影响不同.仿真算例验证了模型和算法的有效性.     

超低轨道地球卫星脉冲星/星光折射/光谱组合导航

针对超低轨道地球卫星导航自主需求,提出了一种脉冲星/星光折射/光谱测速组合天文导航方法。首先根据地球超低轨道卫星运行轨道动力学方程建立导航系统状态模型;分别根据脉冲到达时间差和星光折射角与天体光谱频率建立导航系统量测模型;使用Unscented卡尔曼滤波方法,降低随机误差对导航精度的影响,使用基于UKF的信息融合方法,有效融合了三种天文导航方法结果数据。经计算机仿真分析,该组合导航方法位置导航误差均值为85.62m,速度误差均值0.190m/s,能够满足超低轨道地球卫星在轨运行导航需求。     

自由返回轨道的P平面参数快速设计方法

针对自由返回轨道求解过程中地心轨道类型变化造成的B平面参数方法计算失败问题,提出一种基于P平面参数的自由返回轨道快速设计方法。首先,基于轨道半通径参数的普适性,给出了不同轨道类型的P平面参数定义,建立了以P平面参数为求解目标量的自由返回轨道求解模型。其次,给出了基于P平面参数的自由返回轨道快速设计方法,在构建的瞬时地月惯性系下,以平面双二体自由返回轨道作为初值,实现了高精度力模型下的自由返回轨道快速求解。对8种构型自由返回轨道的设计结果表明,P平面参数具有类似于B平面参数的大收敛域,且有效解决了轨道类型变化对计算的影响,可直接应用于中国后续月球探测任务轨道设计。     

改进PointNetLK的点云智能配准与位姿图优化方法

针对空间在轨服务任务中的非合作目标相对位姿测量问题,提出一种目标可测部位点云的智能配准方法。首先,通过Straight Through滤波算法对半物理仿真平台采集得到的点云进行目标提取,以消除背景数据等杂乱信息;其次,改进PointNetLK神经网络点云配准算法,将提取后的点云数据作为输入,从而获得初步配准结果,解决非合作目标先验信息缺失导致的无法配准问题;最后,建立基于位姿图的优化模型,以降低配准误差,提高配准精度。实验结果表明,与传统迭代最近点(ICP)算法相比,配准综合误差从6.3598降低到1.7291,精度提高约 72.81% 单次耗时从33.16 s降低到4.2 s,效率提升约87.33%,与当前SM ICP等其他算法相比,也具有一定的优势。     

中国深空网首次△DOR联合测轨试验分析

通过分析中国深空网首次△DOR（Delta Differential One way Ranging,双差分单向测距）联合测轨试验的干涉测量事后数据,重点从观测量随机精度、闭合时延等方面讨论了国内深空网与国内VLBI（Very Long Baseline Interferometry,甚长基线干涉测量）观测网、国内深空网与国际深空网的联合干涉处理情况,并与ESOC（European Space Operation Center,欧洲空间操作中心）数据处理结果进行了比对.试验结果表明：我国深空网已具备独立或联合开展深空探测器导航测轨的系统支持能力;深空站系统具备高速率数据接收、采集、记录、传输能力,采集数据处理精度优于1 ns;深空网干涉测量信号处理中心具备多体制信号的干涉处理分析能力,其分析精度与ESOC处理精度差异在0.1 ns量级.     

脉冲雷达凝视模式下低轨碎片相对运动特征分析

脉冲雷达凝视模式是一种探测低轨小碎片的重要方式,其中低轨碎片与雷达的径向运动关系是碎片统计分析的基础.为推算碎片相对雷达的运动特征,基于坐标系的相互转换,先利用雷达测站信息,得到可见碎片在地惯坐标系下的轨道信息;再逆向将目标不同时刻的位置转换至测站坐标系下,得到目标的径向运动特征.在不同要素（雷达波束指向、雷达仰角、测站纬度、目标距离）下,仿真及对比分析了各个参数对速度模糊度、雷达可见性等的影响,得出在凝视模式下目标相对雷达近似做径向匀加速运动,这为后续试验工作提供有力的依据.     

磁强计在航天器上的应用与前景

随着科学的进步,磁强计已被广泛地应用于航天器.本文首先根据磁强计测量原理的不同,对其进行分类.分别介绍各类磁强计的物理测量原理,描述其特性、精度、适用范围.同时概括目前在航天器得到较广泛应用的磁强计.在此基础上,进一步具体分析磁强计作为卫星载荷、姿态测量和控制以及自主导航轨道计算的方法、作用和特点.然后,针对这三个方面应用指出其在航天器上应用存在主要问题和关键技术.最后,对磁强计在航天器上的应用进行总结.同时对其未来的发展进行展望,磁强计在航天器上仍有着良好发展前景.     

交会对接寻的段水平双脉冲交会轨道精确求解

基于考虑摄动影响的精确轨道动力学模型,对交会对接寻的段水平双脉冲交会轨道的精确求解方法进行了研究。提出了将控制脉冲的俯仰角近似转化为控制时刻的轨道幅角,从而调整脉冲控制时刻以消除径向速度增量的方法,精确求解首末水平双脉冲的启控时刻;引入导引终点位置偏差的比例控制方法,精确求解水平双脉冲的精确控制量。仿真结果表明,2轮整体迭代可获得首末水平脉冲控制时刻和控制量的精确值,脉冲水平特性达到俯仰角小于1°,导引终点相对位置精度达到10m,验证了该方法的正确性。     

一类边传递图

令G=（α,β｜α^n=β^2=1,α^β=α^r）,r＆^2≡1（modn），是图Г的一个自同构群。目的是研究关于G一边传递图的性质．运用置换群和代数图论的相关理论，获得了这类图的完全分类，它们是一些互不相交的圈和完全二部图的并。     

DORIS geodesy: A dynamic determination of geocentre location

Geoscience Australia contributed a multi-satellite, multi-year weekly time series to the International DORIS Service combined submission for the construction of International Terrestrial Reference Frame 2008 (ITRF2008). This contributing solution was extended to a study of the capability of DORIS to dynamically estimate the variation in the geocentre location. Two solutions, comprising different constraint configurations of the tracking network, were undertaken. The respective DORIS satellite orbit solutions (SPOT-2, SPOT-4, SPOT-5 and Envisat) were verified and validated by comparison with those produced at the Goddard Space Flight Center (GSFC), DORIS Analysis Centre, for computational consistency and standards. In addition, in the case of Envisat, the trajectories from the GA determined SLR and DORIS orbits were compared. The results for weekly dynamic geocentre estimates from the two constraint configurations were benchmarked against the geometric geocentre estimates from the IDS-2 combined solution. This established that DORIS is capable of determining the dynamic geocentre variation by estimating the degree one spherical harmonic coefficients of the Earth’s gravity potential. It was established that constrained configurations produced similar results for the geocentre location and consequently similar annual amplitudes. For the minimally constrained configuration Greenbelt–Kitab, the mean of the uncertainties of the geocentre location were 2.3, 2.3 and 7.6 mm and RMS of the mean uncertainties were 1.9, 1.2 and 3.5 mm for the X, Y and Z components, respectively. For GA_IDS-2_Datum constrained configuration, the mean of the uncertainties of the geocentre location were 1.7, 1.7 and 6.2 mm and RMS of the mean uncertainties were 0.9, 0.7 and 2.9 mm for the X, Y and Z components, respectively. The mean of the differences of the two DORIS dynamic geocentre solutions with respect to the IDS-2 combination were 1.6, 4.0 and 5.1 mm with an RMS of the mean 21.2, 14.0 and 31.5 mm for the Greenbelt–Kitab configuration and 4.1, 3.9 and 4.3 mm with an RMS 8.1, 9.0 and 28.6 mm for the GA_IDS-2_Datum constraint configuration. The annual amplitudes for each component were estimated to be 5.3, 10.8 and 11.0 mm for the Greenbelt–Kitab configuration and 5.3, 9.3 and 9.4 mm for the GA_IDS-2_Datum constraint configuration. The two DORIS determined dynamic geocentre solutions were compared to the SLR determined dynamic solution (which was determined from the same process of the GA contribution to the ITRF2008 ILRS combination) gave mean differences of 3.3, −4.7 and 2.5 mm with an RMS of 20.7, 17.5 and 28.0 mm for the X, Y and Z components, respectively for the Greenbelt–Kitab configuration and 1.1, −5.4 and 4.4 mm with an RMS of 9.7, 13.3 and 24.9 mm for the GA_IDS-2_Datum configuration. The larger variability is reflected in the respective amplitudes. As a comparison, the annual amplitudes of the SLR determined dynamic geocentre are 0.9, 1.0 and 6.8 mm in the X, Y and Z components. The results from this study indicate that there is potential to achieve precise dynamically determined geocentre from DORIS.